Can surface cracks and unipolar arcs explain breakdown and gradient limits? 
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We argue that the physics of unipolar arcs and surface cracks can help understand rf breakdown, 
and vacuum arc data. We outline a model of the basic mechanisms involved in breakdown and ex- 
plore how the physics of unipolar arcs and cracks can simplify the picture of breakdown and gradient 
limits in accelerators, tokamaks as well as laser ablation, micrometeorites and other applications. 
Cracks are commonly seen in SEM images of arc damage and they are produced as the liquid metal 
cools, they produce the required field enhancements to explain field emission data data and can pro- 
duce fractures that would trigger breakdown events. Unipolar arcs can produce currents sufficient 
to short out rf structures, should cause the sort of damage seen in SEM images, should be unstable 
and possibly self-quenching as seen in optical fluctuations and surface damage. 

PACS numbers: 29.20.-c, 52.80.Vp 



I. INTRODUCTION 

In this paper we explore how well arcing can be ex- 
plained by the properties of surface cracks and unipolar 
arcs. The question of explaining breakdown, arcing and 
gradient limits presents a unique problem, since there are 
over 100 years of reliable published data on an enormous 
variety of phenomena that seem related, but no simple 
explanation has been adopted that can easily be applied 
to clarify or predict the overall physics [TH5]- While it 
is always possible to introduce a variety specific mech- 
anisms that can be narrowly applied, these may not be 
generally useful to explain or predict specific results. We 
find that the mechanisms of unipolar arc physics, com- 
bined with surface cracking, can explain a significant frac- 
tion of the data, and a further analysis seems useful, how- 
ever these mechanisms are not mentioned in most of the 
literature on arcing. Although most of our examples are 
from rf breakdown, the conclusions should have wider 
applicability. The ultimate test of a model is whether 
the ideas are simple, complete and general enough to be 
useful. This paper is an outline of these ideas. 

Our picture of arcs is summarized in Fig. 1 6, 7J. We 
argue that two processes seem to control arcing: 1) the 
formation and fracture of cracks and small structures, 
and, 2) the evolution and properties of unipolar arcs . 
Theoretically, we divide the arcing process itself into four 
elements: 1) mechanical failure of the surface, producing 
fragments, 2) initial ionization of the fragments by field 
emission (FE) currents, 3) exponential density growth 
of the unipolar arc to some equilibrium state, and, 4) 
surface damage produced by the arc. The unipolar arcs 
act as virtual cathodes and produce currents that short 
the rf cavity or other high gradient structure. 

The study of these phenomena has been complicated 
by the speed and unpredictability of the arcs, as well as 
the large dynamic range of the experimental parameters 
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FIG. 1: The arc process is controlled by: 1) fracture of 
high field areas at crack junctions and, 2) the evolution of the 
unipolar arc driven by sheath parameters. 



and the numerical complications involved in simulations, 
where the densities involved seem to exceed the applica- 
bility of the Particle-in-Cell (PIC) codes used for most 
plasma calculations. The difficulties involved in accu- 
rately modeling plasma / surface interactions for very 
dense plasmas are a significant limitation on modeling 

i- 

Numerical modeling of the initiation of the arc using 
a Particle-in-Cell (PIC) has been described in a number 
of papers [3 [7]. Once an arc starts, the surface electric 
field and field emission increase, increasing ionization of 
neutrals, causing an increase in the plasma density. This 
density increase decreases the Debye length and causes an 
increase in the surface electric field, ultimately produc- 
ing an exponential increase in both the electric field and 
density, with time. PIC simulations of the unipolar arc 
model for vacuum arcs relevant to rf cavity breakdown 
show that the density of plasma formed above the field 
emitting asperities can be as high as 10 26 m -3 . The tem- 
perature of such plasma is low, in the range of 1 — 10 eV. 
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While we find that the basic mechanisms can be de- 
scribed simply and some results can be evaluated easily, 
however obtaining more precise results using numerical 
modeling is complicated by the multidisciplinary nature 
of the problems. We will describe some of the basic mech- 
anisms, simple results and more difficult calculations. We 
show how these arguments apply to questions about field 
emission, breakdown, nonideal plasmas, plasma instabil- 
ities and quenching, arc suppression, gradient limits, fre- 
quency dependence, magnetic field effects, etc.. 

II. ELEMENTS 

Since the literature on both surface cracking and 
unipolar arcs in this context is somewhat limited we re- 
view the the relevant physics of these phenomena. 

A. Arc evolution 

We have described how arc evolution can take place 
in rf structures and other environments, see Fig lb. We 
assume that the overall process can be divided into four 
more or less independent stages, surface failure, initial 
ionization, plasma evolution and surface damage. If a 
dense plasma is created on the surface, by laser ablation 
for example, the further evolution of this plasma should 
be expected to be be similar to that of an rf plasma. Sur- 
face failure due to Maxwell stresses has been modeled by 
means of Molecular Dynamics (MD) and the initial ion- 
ization stage and later plasma evolution by means of a 
PIC code. The validity of the results of the PIC code 
is limited in the case of high density, non-ideal plasmas, 
but calculations of nonideal sheath plasmas have been 
done with MD codes. We assume the limiting gradient 
for any system will be determined by a combination of 
the surface damage, which determines the local field en- 
hancements, 

Elocal — P E avera g e 

and the surface failure mechanism. 



B. Unipolar Arcs 

Unipolar arcs were first described by Robson and 
Thonemann in 1959 as an explanation for the existence 
of isolated cathode spots on metal surfaces immersed in 
the plasma of a gas discharge |8J . Unipolar arc phenom- 
ena received extensive study and analysis in the 1970's 
and 1980's as the primary mechanism that determined 
the impurity content of limiter tokamaks. Schwirzke and 
others described both experimental and theoretical work 
with these arcs [9j. As more tokamaks were built with 
divertors, however this mechanism seemed to become less 
relevant, although that may be changing as the physics of 



the ITER tokamak is better understood [10] . The most 
recent and thorough study of unipolar arcs is being done 
by Kajita, who uses laser ablation to produce a plasma 
on a metallic surface that starts the unipolar arc phe- 
nomenon |llj . 

Unipolar arcs can than travel freely on the surface or 
be guided by a magnetic field in the characteristic ret- 
rograde motion that has been identified in many exper- 
iments. The high plasma densities are associated with 
a large plasma pressure which should be responsible for 
particulate production. We have found that the interface 
between the plasma and the liquid surface can become 
turbulent due due to the high plasma pressure in a dense 
arc and the scale of the turbulence is a function of the 
plasma pressure. Some of the parameters of the unipo- 
lar arc plasma can be experimentally estimated from the 
dimensions of the damage produced and measurements 
from SEM images imply the density is very high. 



C. Cracks 

Arrays of cracks are seen in many SEM images of arc 
damage. We believe these cracks are the result of the 
cooling of the melted surface that takes place in two 
stages; first cooling from high temperatures to the so- 
lidification point of the metal, followed by cooling from 
the melting point to room temperature, where the solid 
contracts by an amount Ax — xaAT, where T is the 
temperature, x represents the dimensions of the damage 
and a is the coefficient of linear expansion. 




FIG. 2: SEM image of the center of an arc damage crater. 
The image shows both cracks, with many crack junctions, and 
smooth structures characteristic of a chaotic surface smoothed 
over by surface tension. There is a wide variety of structures 
seen in arc damage, and this image is selected to show both 
cracks and evidence of turbulent structures that have been 
smoothed by surface tension. This image should not be con- 
sidered "typical". 
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During the liquid cooling phase, surface tension would 
smooth the surface, and the relation between the cooling 
time and the scale of surface irregularities seen in SEM 
images can be estimated from the dispersion relation, 

uj 2 = a\k\ 3 /p 

where u>, a and p are the frequency, surface tension con- 
stant and density of the liquid metal, and k is the wave 
number [12, 13]. For copper structures, where we assume 
that thin heated volumes sit on essentially cold surfaces, 
the thermal contraction is approximately 2% of the di- 
mensions of the melted area. We find that the typical 
cooling time constants are in the range of a few hun- 
dred ns for accelerator cavities and the structures seen in 
SEM images of rf cavity damage have radial dimensions 
on the order of a few microns. The two stage cooling 
process seems to result in SEM surfaces that are some- 
what smooth at the 1 micron level, but contain cracks 
with sharp edges at the 1-10 nm level that cover ^2% 
of any large solidified area of copper. 

liquid response time 
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FIG. 3: The relation between the cooling time and structure 
radius for liquid metals compared with data from 805 MHz 
cavity arcs and unipolar arcs in coaxial lines. 



III. DETAILS 
A. Field Emission 

The process of vacuum breakdown was identified by 
students of Michaelson and Millikan almost 110 years 
ago and field emission was of the first mechanisms to 
be described using quantum mechanics by Fowler and 
Nordheim in the 1920's when it was discovered that field 
emission currents were proportional to the applied elec- 
tric field raised to a high power, i ~ E n , with n around 
14 at high surface fields [HLH]. We find that the standard 
method of analysis using Fowler-Nordheim (FN) plots to 



estimate the field enhancement factor can be unnecessar- 
ily abstract, and yields a number (the enhancement fac- 
tor j3) that has little fundamental importance. We prefer 
to plot both the experimental data (currents, radiation 
levels, etc.) against electric field on a log-log plot, along 
with the FN predictions, although the space charge limit 
and other experimental parameters can also be displayed. 
Using this method, the two lines are offset by factors that 
can measure the total emitter area, duty cycle, enhance- 
ment factor and corrections due to the cavity geometry. 
Because field emission currents depend on the electric 
field raised to a very high power, and enhancement fac- 
tors are somewhat difficult to measure experimentally, 
few measurements of the area of field emitters are in the 
literature. 
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FIG. 4: Field emission can be plotted to display the vari- 
ables associated with its measurement. The horizontal and 
vertical offset of the data and theory curves are essentially 
the enhancement factor and emitter area but the effects of 
thermal emission, different work functions, duty cycle, struc- 
ture geometry etc., as well as systematic and statistical errors 
in measurement can also be displayed graphically. 



Following early work by Dyke et. al. showing that 
Ohmic heating of tungsten needles could produce break- 
down, combined with considerable evidence that the re- 
quired field enhancements and current densities could be 
produced with cylindrical asperities with rounded ends, 
Ohmic heating was widely accepted as an explanation 
for breakdown, although asperities of the expected di- 
mensions were not found [? ]. We have shown that 
cracks, more specifically crack junctions, can provide the 
required field enhancements, and emitter areas (when 
many of them are added together) to explain field emis- 
sion data. 

Field emission measurements in rf cavities and high 
gradient structures have been made and reported in a 
number of references |16j . These measurements assume 
that the number of emitters is known and understood. 



4 



We assume that a large number of much smaller, localized 
emitters at crack junctions, contribute as one emitter. 

B. Enhancement Factors 

Following Feynman, we describe the surface field of a 
conductor as a function of the local curvature of the sur- 
face, comparing the fields at any two points a and b will 
give the relation E a /Eb = rf,/r a , where r is the three 
dimensional radius [T7]. Small radii give high fields. We 
find these small radii at crack junctions, where the radii 
are too small to be resolved by SEM optics. Numerical 
analysis has shown that these crack junctions can pro- 
duce enhancement factors in the range of f3 ~ 200. 



* le-6 




FIG. 5: The crack junction mesh. 

In another example of field enhancements, one can de- 
scribe the sheath potential of a tenuous plasma as an en- 
hancement of an applied field that would add the field in 
the sheath to the externally applied field. As the plasma 
density increases, this picture becomes less appropriate. 

C. Breakdown Mechanisms 

It has been shown that breakdown occurs at fields 
near 10 GV/m. These local fields would cause mehanical 
stresses on the order of, 

a = e Q E 2 /2 = AA x 10 s MPa, 

which would be pulsed and subject to fa- 
tigue, accompanied by high field emission current 
densities, and perhaps local heating. Under these condi- 
tions mechanical failure would be expecte. Electrostatic 
fracture, Ohmic heating, electromigration, fatigue and 
creep can all explain the the mechanical failures that 



could trigger breakdown, for example the breakdown 
rate BDR ~ E 30 behavior seen in some experiments. 
For example, both Ohmic heating and electromigration 
should be proportional to the current density squared, 
since field emission produces current densities in the 
range j ~ E 14 . Electrostatic fracture is similar to 
field evaporation, which is governed by processes that 
produce field scaling from rate ~ 2? 30-150 , j-j^ seems 
compatible, although the effects of creep and fatigue are 
not well understood. The more difficult problem is to 
find a mechanism that is compatible with the damage 
seen in SEM images. 

Alltlhough we favor the model of electrical stress and 
fatigue as a trigger for breakdown, we find that the BD 
mechanism itself is less interesting than a description of 
the environments that are highly stressed in many dif- 
ferent parameters. It is important to understand nature 
of asperities and their geometry to understand if it is 
possible to suppress breakdown. 



D. Parameters of Nonideal Plasmas 

Simulations with PIC codes have shown that field emis- 
sion, significant sheath potentials, high densities of neu- 
trals, along with self-sputtering can produce an environ- 
ment where the density rises essentially exponentially 
while the electron and ion temperatures remain relatively 
low. This increasing density is associated with a decreas- 
ing Debye length, 

\ D = \J eok B T/n e e 2 , 

so the number of particles in the Debye sphere eventually 
becomes less than one and the nonlinearity parameter, 
, which measures the ratio of the electrostatic poten- 
tial energy divided by the kinetic energy of the plasma, 
becomes large [18]. 

Recent numerical analysis of high density, non-ideal 
plasma sheaths has shown that for high density plas- 
mas the properties of the plasma can be estimated us- 
ing molecular dynamics. The corrections to simple esti- 
mates of sheath potential, Debye length and surface elec- 
tric field required by the non-ideality condition due to 
the high densities involved are not large. Plasma densi- 
ties were estimated from the scale of damage, where tur- 
bulence of produced by the plasma pressure is balanced 
against the smoothing produced by the surface tension. 
Since it is difficult to know the cooling time with much 
precision, these measurements function as an upper limit 
on the scale of turbulence and a lower limit on the plasma 
density. This procedure produces estimates of the sur- 
face plasma density n ~ 10 25 m -3 , surface electric fields 
E ~ 2 x 10 9 GV/m, for electron temperatures of 10 eV 
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E. Space Charge Oscillations 

PIC codes have shown that when held emitters can 
ionize dense gas near the surface, a positively charged 
plasma is produced, and the sheath potential of the 
plasma that is created increases the field on the field 
emitters until they become space charge limited. 

The space charge limit for continuous currents between 
two plates is expressed using the Child-Langmuir Law, 

4e p^SV^ 

where / is the anode current, the current density, and 
S the anode surface inner area [19] . While the Child- 
Langmuir Law applies to thermionic emission, the appli- 
cation of this idea to field emission is not entirely straight- 
forward. Thermionic emission of electrons is essentially 
constant, with fluctuations governed by variations in the 
temperature of the emitter. With field emission, how- 
ever, the current density is proportional to the electric 
field to some high power, (i ~ E 14 ), thus, fluctuations 
in the electric field will instantly alter the emitted cur- 
rent density, and fluctuations in the density of emitted 
electrons will immediately alter the electric field. These 
processes can produce fluctuations. 

PIC code results show that the space charge limited 
current is not continuous on a microscopic scale. Elec- 
trons are emitted from the surface in bunches, they move 
a few microns away from the cathode where they produce 
a negatively charged electron cloud that erodes due to 
electrons moving both toward and away from the cath- 
ode. This repetitive behavior produces an oscillation in 
the field emitted current at a frequency of about 1 THz 
[20] . We are not aware of experimental observation of 
this phenomenon. 



F. Plasma Fluctuations and Quenching 

Vacuum arcs can be unstable. Fluctuations in the op- 
tical emission of arcs have been recorded in streak camera 
experiments and one of the defining properties of unipo- 
lar arcs is their random, discontinuous, trail of surface 
damage. Optical fluctuations occur at frequencies up to 
a few hundred MHz. We describe the fluctuations seen in 
unipolar arcs as a similar mechanism to the fluctuations 
in the space charge limited field emission described above 
in Section IIIA. 

Although over long time scales the plasma should 
maintain quasi-neutrality, the mechanisms controlling 
the electron and ion densities are quite different, and 
have different timescales. Since the field emitted cur- 
rent density will be proportional to, %fe ~ E le , the field 
emission current will respond instantly and nonlinearly 
to changes in the surface field, and the electrons can be 
rapidly thermalized in a dense plasma. 



We assume that the fundamental ion density increase 
is governed by self-sustained self-sputtering, 

a/37 > 1) 

where a is the probability that a sputtered cathode atom 
becomes ionized, j3 is the probability that the ionized 
atom returns to the cathode, and 7 is the sputtering yield 
[2]. The ion density, rii, should respond slowly, since 
the time constant for density changes would depend on 
collisional diffusion [21] in the arc, 

dnjdt = D V 2 n 2 , 

which is a function of the density, rii , since the diffusion 
constant is inversely proportional to the plasma density, 

D = vth/iv — 1/m. 

As the arc evolves and the density increases, the large, 
dense arcs should become more stable to ion density fluc- 
tuations, with time constants proportional to, Tj ~ rii, 
the time constant for field emission, however, should not 
change and the electron thermalization time should be- 
come shorter as the density increases like, r e ~ 1/rt;. 
This difference between the ion and electron density sta- 
bility could complicate the ability of the plasma to main- 
tain quasi-neutrality under rapid high current field emis- 
sion. 

As the arc evolves, surface fields created by the sheath 
potential become large enough to produce field emission 
currents that can short out the sheath potential and lo- 
cally quench the arc before quasi-neutrality can be estab- 
lished. As shown in reference [18] . the current required 
to short out a plasma sheath in time At is equal to, 

i s = e E/At, 

for times of At = 1 ns, the required current would be ~30 
MA/m, which could be produced by a field of ~3 GV/m, 
which is compatible with simulations produced by both 
PIC and MD codes. The remaining dense plasma is then 
either able to restart the arc nearby, or, after a time re- 
quired to equilibrate the locally dense plasma, restart in 
the same location. These densities are compatible with 
data taken with 805 MHz rf structures, which have arc 
damage diameters of 0.5 mm diameter and shorting cur- 
rents on the order of 10 A. This argument seems to pre- 
clude plasma / surface fields significantly larger than 3 
GV/m and current densities larger than 30 MV/m. 

The comparatively low current density of 30 MA/m 2 
is not large enough to produce significant Ohmic heating 
of the surface. This seems to conflict with the current 
densities required by the ecton model of G. Mesyats [5]. 
In that model, current densities of ~ 10 13 A/m are re- 
quired to produce a local Ohmic heating explosion of the 
liquid metal that maintains the arc. 
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G. Frequency Dependence of Gradient Limits 

The maximum operating gradient of a given rf or DC 
system could operate should be a function of two vari- 
ables; 1) the maximum local field at which the surface 
would fail, due to tensile stresses, heating, electromigra- 
tion, fatigue or some other effect, and, 2) the overall de- 
sign of the system itself, which determines the stored 
energy deposited through the arc, the way power is ap- 
plied, discharge length, the way the power to the arc is 
turned off (suddenly or slowly) which all seem capable 
of affecting the surface damage, and ultimately the field 
enhancements seen by the surface [22] . 

A large body of data showed very early that DC break- 
down occurred when local fields reached 7-10 GV/m 
over many orders of magnitude variations in the gap 
length. Although there are not many rf measurements, 
data also show that this threshold also seems to apply to 
systems around 1 GHz. These results imply that there 
is no frequency dependence to high gradient breakdown 
as a function of the local electric field, Ei oca i. On the 
other hand, there is an extensive literature that show 
that higher frequencies achieve higher gradients. 

H. Magnetic Fields 

We have found that under some conditions the maxi- 
mum rf electric field that can be maintained in the pres- 
ence of an externally applied electric field is reduced from 
that seen without the external magnetic field. 

The beam optics of field emitted beams in magnetic 
fields have been studied experimentally. The beams were 
actually found to be hollow, with a radius that was di- 
rectly proportional to the applied electric field and in- 
versely proportional to the square of the static magnetic 
field. This is consistent with a picture of field emission 
from a foil-less diode, where the electric and magnetic 
fields were not parallel. 

SEM images of arc damage in copper with a magnetic 
field shows that surface cracking is confined to a very 
small area (10 - 20) /im in diameter surrounded by a much 
larger area that shows signs of being melted. We assume 
that the radial growth of the arcs were confined by the 
magnetic field, and when the arcs cooled they cooled from 
the outside in, leaving the last metal to solidify to absorb 
all the thermal contraction. These central damaged areas 
have a much higher crack density, that we associate with 
the lower electric fields that could be maintained on the 
surface. 



I. Other Environments 

Although we consider arcing primarily in the context of 
rf linacs, It is useful to see how generally these arguments 
apply to other environments, such as tokamak first walls 
[TU] . laser ablation [TT], micrometeorite impacts [53] and 



possibly such examples as electron beam welding. Arc 
damage in laser ablation, tokamak first walls and mi- 
crometeorite impacts seems essentially identical. 

The study of unipolar arcs is not an active field. Al- 
though vacuum arcing phenomena are seen in many en- 
vironments, there is little contact and little coordination 
between different approaches used in these fields. Arc 
damage seems to be quite similar between micrometeorite 
impacts, laser ablation targets and tokamak first walls, 
however the lack of a model or a common approach to 
understanding the basic mechanisms at work in unipolar 
arcs has further slowed progress. 



IV. USEFUL EXPERIMENTS 

There have been 110 years of experimentation on vac- 
uum arcs, most of them guided, to some extent, by mod- 
eling and theory, nevertheless there is still disagreement 
about the nature of these arcs and the mechanisms that 
drive them. We believe the reason for this situation is 
that the arcs are small and unpredictable, and many pa- 
rameters (which are individually hard to measure) evolve 
very rapidly over a many orders of magnitude. While 
models exist, theory and modeling are complicated by 
the large number of mechanisms that seem to be involved 
in arc evolution and high density plasmas, that require a 
complicated, non-Debye analysis of even basic properties. 

There are a number experimental directions that could 
prove promising. Space charge oscillations have not been 
seen as far as the authors are aware, although they would 
be somewhat difficult to detect because of the high fre- 
quency involved. Likewise, quench of unipolar arcs has 
never been studied in any detail, as it is difficult to access 
the surface underneath the plasma. And lastly, it would 
be useful to have more systematic data on the damage 
mechanism in normal arcs. While this field has produced 
considerable data, the problem is that little of it was sys- 
tematically selected and documented. 



V. CONCLUSIONS 

We have demonstrated that the physics of unipolar 
arcs and surface cracking seems highly relevant to the the 
phenomenon of arcing. Crack junctions can provide the 
high field enhancements seen in experimental studies of 
field emission and breakdown, they are formed naturally 
as arc damage cools and they are capable of triggering 
breakdown events when they fracture. Likewise, unipolar 
arcs are the prototypes of single sided arcs, that can func- 
tion as cathode spots. We find that the physics of these 
objects, which has not received much specific attention 
is relevant to many fields. 

Although unipolar arcs have been called ubiquitous, 
the literature on this phenomenon, both experimental an 
theoretical, is not extensive. We believe one reason for 
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this is that the dense plasmas and plasma / surface in- 
teractions require very specific techniques to cope with 
the nonideal (non-Debye) plasmas, that are not well ad- 
vanced. 

Wc have shown that cracks are commonly seen in SEM 
images of arc damage and described how they are pro- 
duced as the liquid metal cools below the melting point 
to room temperature. We have shown that cracks can 
produce the required field enhancements to explain field 
emission data data and can produce fractures that would 
trigger breakdown events. Although unipolar arcs, and 
non-ideal plasma surface interactions are not well stud- 
ied, we have shown that field emission of electrons pro- 
duced in the plasma sheath can produce currents (1 - 
1000 A) sufficient to short out rf structures, should cause 
the sort of damage seen in SEM images. These plasmas 
should be unstable and possibly self-quenching as seen in 
optical fluctuations and surface damage in a variety of 



experiments. 

Although the internal structure and evolution of arcs 
and arcing has not evolved to any unanimity in the field, 
wc find that the physics of cracks and unipolar arcs seem 
highly relevant and perhaps fundamental to these phe- 
nomena. 
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